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Abstract 

We revisit the problem of semi-flexible chain condensation by neutral 
depleting agents {e.g. colloidal spheres or flexible polymers) by using a 
simple formalism that allows us to address its main features without spec- 
ifying the kind of depleting agents. Correlations between depleting agents 
are shown to produce a reswelling of the chain at high enough volume frac- 
tion, consistent with an earlier analysis by Grosberg et al. {Biopolymers, 
21:2413, 1982) in the context of DNA condensation induced by neutral 
polymers. It is also shown that the conditions under which spherical col- 
loids can condense a semi-flexible chain are strongly restrictive, unlike 
what has been recently claimed. The formalism allows us to compare the 
efficiency of colloids vs polymers as condensing agents. 



1 Introduction 

Since the early and pioneering work of Asakura and Oosawa^ in colloid-polymer 
mixtures, depletion forces have been invoked to explain a number of behaviors 
in various systems. In the original system considered by Asakura and Oosawa, 
for example, or in the case of an asymmetric binary mixture of hard spheres, 
the effective interaction between pairs of large spheres is attractive at short 
distances due to the depletion of the s maller component in between and can 
induce fluid-fluid phase separation.l^Eni The theoretical description of depletion 
potentials has even recently become still more quantitative for the simplest 
geometries like the asymmetric hard sphere mixturcEI In particular, it has been 
realized that accurate depletion potentials must take into account correlations 
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among depleting agents. Those correlations give rise to an energy barrier in 
the effective interaction between pairs of large colloids that they have to cross 
before feeling the classical attractive (Asakura-Oosawa) potential. This feature 
strongly affects the phase diagram of asymmetric mixtures.^ 

Depletion forces are also responsible for the polymer salt induced (psi) or 
ip condensation at the single chain level: each isolated chain of semi-flexible 
polyelectrolyte is condensed into a dense toroid above some critical concen- 
tration of neutral polymer and salt. Since its discovery in the seventies by 
Lerman, in the particular case of DNA and polyethyleneglycol (PEG), this 
condensation has been addressed extensively both experimentally and theo- 
retically.'' ^■^J^'i^*^^^ Its physical origin is attributed to depletion 
interactions between DNA segments induced by PEG; the volume occupied by 
the DNA depletes a certain amount of PEG molecules and therefore induces 
an effective attraction between DNA segments. The theoretical methods used 
so far in the context of single chain condensation are usually taken from poly- 
mer physics theories. Due to the connectivity of chain segments, correlation 
effects are known to be very important both at the single chain and solution 
level.ESl Consequently, the emphasis has mainly been put s o far on these corre- 
lations when describing the chain coUapsei ^'^l^^ l ^^ l -'^^ l ^^ l ^'-'l^-'^ l In particular, it has 
recently been sugges ted that spherical colloids of appropriate sizes could also 
condense DNA.EHUni 

In this spirit, Odijk proposed a model for the nucleoid of 
bacterial cells,'23Sl wherein cellular DNA is condensed by the depletion inter- 
actions between DNA segments induced by proteins of the cytoplasm that bear 
the same charge as DNA. In all these works, the correlations between colloids 
were limited to the second virial level; also, no restrictions on the colloid size 
were necessary to obtain the DNA condensation. However, as first emphazised 
by Grosberg et ai, the omitted correlations between PEGs in the case of ■0-DNA 
condensation can qualitatively change this pictur^l - above another, higher, 
critical concentration of PEG the DNA condensate reswells. One speaks then 
of a 'reentrant condensation'. This phenomenon has indeed been observed by 
Vasilevskaya et al. by fluorescence microscopy,!^ and is very similar to what is 
predicted and observed for gels immersed in polymeric solution where there is a 
collapse and reswelling of the gel for two different critical polymer concentrations 
in the solutionPSl 

Apart from the work of Grosberg et al., the level of description of depletion 
effects in single chain condensation contrasts with that in the recent literature 
on the phase behavior of colloid-polymer mixtures, which latter is mainly de- 
scribed within the formalism of simple liquid theories where the correlations 
between depleting agents are highlighted.!^ As will be shown in the present 
paper, insufficient account of correlations in the single chain collapse problem 
can lead to missing important effects. In particular, inclusion of correlation ef- 
fects between depleting agents in the case of semi-flexible chain condensation by 
neutral polymers leads to the requirement of a minimal degree of polymerisation 
for the polymer to be able to condense the chain, and for re-entrance.— This 
prediction has been checked experimentally. In the case where spherical colloids 
are the depleting agents, as discussed below, the effective interactions between 
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chain segments might not become attractive at all, whatever the concentration 
of colloids used. 

The aim of this work is to revisit the effective interactions between chain 
segments induced by depleting agents, by mixing polymer physics tools and 
simple liquid theories. We propose a simple formalism that accounts for the 
main features of semi-flexible chain condensation and treats at the same time 
a high level of correlations between depleting agents. In the next section, we 
present the basic framework of our model, without specifying the kind of de- 
pleting agents used. Then the application of this model and some quantitative 
results for hard spheres are given in the third section. Finally, those results are 
discussed and compared to previous approaches to semi-flexible chain collapse, 
and in particular DNA condensation. Moreover, the relative efflciency of flexible 
polymer vs colloid for chain condensation is compared. 

2 General results 

2.1 General framework for the description of semi-flexible 
chain condensation 

The semi-flexible chain will be modeled as a freely jointed chain, of segment 
length Ip and width D.^^ Each segment can be divided into n "monomers" of 
size D and the total number of monomers of the chain is N, considered to oc- 
cupy a spherical volume V. We neglect any kind of intramolecular orientational 
ordering due to the anisotropic shape of the segments, an approximation valid 
in the limit of inflnite chain length as long as the density of the condensate does 
not exceed the critical density of the isotropic- nematic ordering transition of un- 
connected segments.!^ Since our goal in this work is not to draw a precise phase 
diagram of chain condensation, but rather to study the effect of correlations on 
the induced depletion interactions, the choice of an isotropic condensate with an 
imposed spherical geometry seems well suited for the tractability of the model. 
In particular, one should be able this way to address the geometric properties 
of depleting agents (degree of polymerization or colloid size) required to achieve 
condensation in a finite range of depletant concentration. 

The partition of the neutral depleting agents inside and outside this volume, 
and their interactions with chain segments, can give rise to effective attractive 
interactions. The volume fraction of depleting agents inside and outside the 
chain volume are respectively yi„ and y. To be more specific with the notation, 
we consider explicitly the case of spherical colloids as condensing agents, but the 
results that we derive in this section will be easily generalized to other depleting 
agents like flexible polymers. The direct interaction (per unit volume) between 
chain segments and colloids is chosen to have the form: 

V Tra-^/D 
where Vcross = f ^p(^ + 'i)^ is the second cross virial coefficient between chain 



3 



segment and depleting agents, cg is the segment density inside the chain volume 
and a is the diameter of the colloids. This choice will be justified when we 
discuss the results of our model in sectional All energies are measured in units 
of kT = 1. Let us introduce the reduced concentration X = VcrossCs- The 
grand potential describing the equilibrium properties of the system is 

n[V, y,„] = V |f(2;,„) + F,nt{y^n) + p{y) ~ ^^^^} + FpoiiX) (2) 

The quantities F{y),p{y) and ^{y) are respectively the free energy density, os- 
motic pressure and chemical potential of the pure colloid system. We do not 
specify right now the form of these quantities, but the Scaled Particle Theory 
(SPT) for hard sphere#^ will be used to evaluate them in the next section. 
Similarly we do not specify yet the free energy, Fpoi{X), associated with the 
pure semi-flexible chain. It will consist of excluded volume effects, and confor- 
mational entropy (see Eg ITT)) . 

The system described by the grand potential f2 is in equilibrium with a 
reservoir of colloids and a reservoir of volume, which means that the colloid 
chemical potential and the pressure of the system are imposed. The equilibrium 
equations for the "inside" and "outside" condensing agents are the equality of 
their chemical potential and the equality of osmotic pressure: 

Kyin)+X = n{y) (3) 

p{yrn)-p{:y) + X^i^ + Ppoi{X) = (4) 
7ra /D 

with the osmotic pressure of the condensing agents given by p{y) = —F{y) + 
y^Ko^jQ' pressure due to the semi- flexible chain alone is denoted Ppoi- By re- 
placing X in the osmotic pressure balance with its expression from the chemical 
potential equilibrium, the osmotic balance condition can be rewritten 

Pe//(X)+Ppo;(X) = (5) 

In the last equation, the effective pressure or interaction P^f / is defined by 



y.n = P^-^ [/^(y) - X] (7) 



Fe// {X) = F{y) - F{y;n) + Pij^ {y^n " V) (6) 



If the functional form of the condensing agent chemical potential is known, the 
effective interaction can be estimated by using Eqs. I5I6I and 171 Indeed is 
the only input required to make this estimate. 

In order to obtain very general information on the effective interactions 
between chain segments, we can compute their properties as a virial expansion in 
X.ESl The n-th coefficient, i.e. the contribution in X^ in the osmotic balance, 
is obtained by taking the n-th derivative of the effective pressure Peff{X), 
evaluated at X = 0. The n = and n = 1 derivatives vanish identically when 
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evaluated at X ~ 0. But the second derivative of Peff, when added to that 
from the "bare" term Ppoi, and evaluated aX X = 0, gives the effective second 
virial coefficient between chain segments 

With this general form for the effective virial coefficient that depends only 
on the chemical potential of the pure condensing agent system, we can study 
the role of correlations between colloids. If the colloid chemical potential is 
represented as a virial expansion as a first step, for example, one has for hard 
sphere depleting agents 

/i ~ logy + 8y+ 15?/^ 

At' ~ -+8 + 30y (9) 

y 

where the exact second and third virial coefficients for hard spheres have been 
used!^ Now if we stop the expansion at the second virial level for the condens- 
ing agents (i.e., keeping only the first two terms in Eq. [Sjfor /i'), the effective 
second virial of chain segments is a monotonic decreasing function of the col- 
loid volume fraction and becomes negative above a critical value of y. In the 
framework of the coil-globule transition theory, a negative second virial coeffi- 
cient means that the chain is collapsed from a swollen to globular state where 
the correlation of segment density fluctuations is much smaller than the globule 
size.l^Sl Accordingly there should be an unconditional collapse of the semi-flexible 
chain by colloids of any size once the critical concentration is reached. However 
if the third virial coefficient of colloids is taken into account, and more gener- 
ally higher order coefficients, the effective second virial coefficient between chain 
segments is a non-monotonic function of the colloid volume fraction. Therefore 
the precise value of the minimal effective second virial coefficient determines if 
there is actually a collapse of the chain, as will be discussed in the next section. 

If this is the case, then the effective second virial is predicted to become 
positive again for higher volume fraction of colloids: the chain reswells. This 
behavior can be interpreted in the following way: the interactions between col- 
loids become more and more important as their volume fraction is increased, 
so that the interactions between chain segments and colloids are only a slight 
perturbation of the colloidal fluid. Therefore the system tends to come back 
to the swollen state of the chain, where the colloid volume fraction inside and 
outside are similar, since colloids do not lose any translational entropy this way. 

Let us illustrate qualitatively these results in the case where flexible poly- 
mers are the neutral condensing agents. For concentrations below the overlap 
concentration for polymers, the chemical potential is simply jj, ^ logjj, where 
c is the monomer concentration and M is the molecular weight of the polymer. 
For semi-dilute concentrations, the chemical potential is ~ Mc^/^, taking 
into account scaling laws of neutral polymers?^ The derivative of the chemical 
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potential therefore reads 

^ — ifc«c* 
c 

n'^M^c^/^ if c»c* (10) 

The effective second virial coefficient is non-monotonic as in the case of strongly 
correlated colloids. We expect also a reentrant collapse of the chain at high con- 
centration of polymers. This reentrant collapse for a single semi-flexible chain 
interacting with a solution of flexible polymers was predicted for the flrst time 
by Grosberg et in the context of ^-DNA condensation and, as mentioned 
in the introductory section, this reswelling of the DNA has been observed exper- 
imentally by Vasilevskaya et al^^ Note that according to these authors, such a 
reentrant collapse could be taken into account through a concentration depen- 
dent Flory parameter for the flexible polymer, which is equivalent to including 
a higher level of correlations in the system (see below). 

Related results concerning the role of correlations on the effective second 
virial coefficient in mixtures of hard spheres and also colloid-polymer mixtures 
have been obtained previously. As mentioned earlier, one of the early models 
of the depletion potential in colloid-polymer mixtures is due to Asakura, Oo- 
sawa and then VriJil^ZI _ 'Q■^xl ti^e interactions between depleting agents (here 
the polymers) are not taken into account within those models, and the resulting 
depletion potential between colloids is purely attractive. On the other hand, 
more accurate treatments of this system, including interactions between deplet- 
ing agents, give rise to oscillatory potentials with both attractive and repulsive 
parts.l^ In particular, there is a growing repulsive energy barrier to the short 
range Asakura- Oosawa type attraction as the concentration of depleting agents 
is increased. Consequently the effective second virial coefficient between colloids 
is non-monotonic,!^ a feature that is also predicted within our simple model. 



2.2 Explicit model for the semi-flexible chain 

The free energy associated with the pure semi-flexible chain system consists of 
an excluded volume term, and a conformational entropy contribution: 

Ppol{X) = Fev + Fconf (H) 

with 

VS 



^cross 



V 



F, 



The excluded volume term is limited for the sake of simplicity to the second virial 
term, although full self-consistency within the approximations made here would 
require, for example, developing a scaling form for the unperturbed polymer 
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free energy. But Eq. E| suffices, as our aim is to concentrate on the role of the 
effective second virial coefficient in driving the collapse and swelling transitions. 
Since the chain segment is approximated by a rod of width D and of length 
Ip, the virial coefficient is given by vs = -^IpD. The conformational entropy is 
composed of the cost for stretching (first term) and compressing (second term) 
an ideal chain. 

To determine the equilibrium conditions for the chain in a solution of deplet- 
ing agent (volume fraction y), we need to minimize the grand potential fl - given 
by EqEl- with respect to: (i) the number of colloids "inside" the chain volume; 
and (ii) the volume occupied by the chain. The first minimization leads, as 
before, to Eq|3| or - equivalently - [7| The second gives again Eq[Sl but now, 
with EqE]for Fpoi{X) allowing us to evaluate explicitly Ppoi = —{dFpoi/dV), 
keeping all numbers of particles constant. We find 

(12) 

The terms adding to the colloid-induced contributions Peff{X) are: the ex- 
cluded volume associated with the second-virial interaction between chain seg- 
ments (each with volume vs and number density cs = X/vcross)', and the chain 
conformational terms going as d{R^)/dV ~ X^^^ and d{R^^)/dV ~ 

y~^/^ ~ X^/^. As emphasized in the previous section, the induced depletion 
interactions are given by Pg// = piVin) — p{y) + ^Ji"/Q in the osmotic balance, 
with yin{X) determined from Eq[7| 

Within a virial expansion, the general form of the osmotic balance reads 

-AX^'^ + 5X^/3 + {C- C,oii)X^ + D,ouX^ + ... = (13) 

with ^^coW denoting the contribution from colloid-induced depletion effects. 
The solutions of this equation are extrema of the grand potential EqI21 One has 
therefore to be careful with the stability of the solutions of interest. This form 
for the osmotic balance can be used to discuss qualitatively the various features 
of chain condensation and to guide the numerical analysis. The first step of the 
numerical analysis is to invert numerically the function "chemical potential" - 
see Eq[7| so that yin{X) can be used directly in the effective interaction Eall2[ 
without any virial expansion. 

If there are no colloids, the osmotic balance reads —AX^/'^ + BX^^'^ + CX^ = 
0. Under good solvent conditions for the semi-flexible chain, the term in A 
balances the term in C, while in a ^-solvent, the term in A balances the term in 
B. One can easily check that this equation is just the Flory model for the end-to- 
end radius of the chainPSI In the presence of colloids, the depletion interactions 
induced by the hard spheres change the effective interactions among segments. 
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3 Quantitative results for hard sphere-like de- 
pleting agents 

In this section, we consider explicitly the case of hard spheres as depleting 
agents. The only input needed to compute the depletion interaction Peff be- 
tween chain segments in the previous general model are the properties of the 
pure hard sphere system through the chemical potential. Since we are dealing 
with a hard sphere fluid, we can use the thermodynamic quantities derived from 
the Scaled Particle Theory (SPT) which are known to be accurate as long as 
the colloid volume fraction is not too close to the close packing limit ^ The free 
energy density, osmotic pressure and chemical potential within the SPT read 

ny) - -^(:,7T^+^og-^) (14) 



7ra3/6 V2(l l-y 

M?/) l°ST3^+ 2(1 -y)3 (16) 

It is well known from hard sphere theory that the second and third virial com- 
puted through the SPT are exact, while up to the fifth order the differences are 
small (less than 5%). 

Using the analysis of the previous section, the nature (repulsive or attractive) 
of the effective interaction is determined through the effective second virial 
coefficient between chain segments. The result, calculated analytically according 
to EqlHI is shown for example in figureQlfor the particular value a/D = 15 chosen 
by Sear to study DNA condensation by spherical colloids 

As mentioned before, if the pure colloid system is described only up to 
the second virial level, as in the analysis of Sear, the effective second virial 
coefficient between chain segments becomes negative above a critical colloid 
volume fraction (~ 0.17 in this instance); under these conditions, the chain 
(DNA) should collapse. However if the pure colloid system is described up to the 
third order, the minimum of the effective second virial coefficient is positive and 
there is no collapse of the chain. This result is confirmed when the pure colloid 
system is described through the SPT, the minimum of the effective second virial 
coefficient being slightly higher and still positive. Note that the effect of the 
correlations is relevant for moderate values of colloid volume fraction (y « .2) 
so that the description of the hard sphere fluid through the SPT should be 
valid. Also, the volume fraction at which the effective second virial coefficient is 
minimum does not depend on the ratio a/D. It is determined solely by the level 
of correlations between depleting agents, see, e.g. third virial order vs SPT in 
figure n 

It is straightforward to compute the minimum of the effective second virial 
coefficient within the SPT description of the colloids as a function of the ratio 
a/D. The result is shown in figure |21 
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Figure 1: Effective second virial coefficient as function of y. From bottom to top: 
second virial approximation, third virial approximation and SPT approximation 
for the colloids. a/D = 15 



We see that the colloids can induce a negative second virial, only if a/D < 
0.63 or a/D > 24.17. Now for DNA, we have Ip/D ~ 25, so that according to 
this model colloids whose diameters are bigger than the persistence length of 
DNA should condense DNA. Unfortunately, we also expect the model to break 
down in this limit because the form we chose for the cross virial Vcross should no 
longer be appropriate. On the other hand, the range a/D < 0.63 corresponds 
to molecular sizes rather than colloids, so that we do not expect any realistic 
colloids that behave like hard spheres to condense DNA. Nevertheless, as is 
shown immediately below, we still find chain condensation by spherical colloids, 
but under different conditions. 

Note that the relative magnitude of the depletion interactions, i.e., the ratio 
between the induced second virial coefficient and the bare one, scales like 



vsa^fi'{y) 

It is possible to increase the effect of depletion interactions induced by hard 
spheres by decreasing, say, vs. This corresponds to decreasing the solvent 
quality for the semi-flexible ch ain by adding, for example, a certain amount 
of methanol to aqueous solutiorP2l22I in the case of DNA. Of course, this as- 
sumes that the colloids still behave purely like hard spheres, so that the solvent 
does not induce attractive Van der Waals interactions between colloids. De- 
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Figure 2: Mimimiim of the effective second virial coefficient as function of a/D. 



creasing slightly the quality of the solvent such that the isolated chain is at 
least in a theta solvent should change the balance between bare and induced 
interactions. In our model, the solvent quality is taken into account by the 
substitution vs tvs, where the solvent strength r can be varied between 
(theta solvent) and 1 (athermal solvent). Negative values of r correspond to 
a poor solvent, where the question of chain condensation by depleting agent is 
irrelevant. The numerical solution of the equilibrium equations in the limit of 
infinite chain length {N oo) for a/D — 15 and r — 0.67, is shown in figure 
El It is clear from these results that in this case the chain condensation is con- 
tinuous, while the reswelling is discontinuous. Note that the density of chain 
segments in the coil state in the limit of infinite chain length is X — 0. 

On the contrary, if the chain has a finite length, the conformational entropy 
term can balance the depletion interactions. In this case, we expect a discontin- 
uous chain condensation; for a given y, the repulsive BX^/^ term is dominant 
at low X, while the attractive depletion term —CcoiiX^ is dominant at higher 
X. Above a critical colloid volume fraction, the chain segment density jumps 
from the coil value to the globule value. As in the infinite chain length limit, 
the chain condensate reswells to the coil value at still higher volume fraction, 
because of the strong correlations between colloids. The numerical solution of 
the equilibrium equations for a particular range of parameters exhibiting this 
specific behavior is shown on figure^ accordingly, the solvent quality, i.e. t 
has been tuned. Note that the density of the chain segments in the coil state is 
now finite, although its value is very small (cf figure |31 w 0} . 
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Figure 3: Reduced equilibrium density of chain segments X as a function of 
the colloid volume fraction y in the limit of infinite chain length and reduced 
solvent quality. a/D = 15 , Ip/D = 25 , r = 0.67. 



In both cases of finite or infinite length of semi-flexible chain, we find that 
the reswelling to the coil state is discontinuous. This can be easily interpreted in 
terms of effective virial coefficients; above a certain volume fraction of colloids 
higher than the critical value for chain collapse, the effective second virial coef- 
ficient increases with y and becomes positive again while the third virial is still 
negative. Inspection of higher order induced virial coefficients shows that as the 
order of the virial coefficient increases, it changes sign more often as a function 
of colloid volume fraction. This indicates that the use of the virial expansion 
for the effective interactions is quite limited to the location of a critical value of 
the colloid volume fraction for chain collapse, and is not well suited to describe 
precisely the properties of the condensate itself. For example, the reswelling of 
the chain is not located at the value corresponding to the second change of sign 
of the effective second virial coefficient. Note that a similar abrupt reswelling is 
predicted for gels immersed in a solution of polymers.!^ This is not surprising 
due to the similarities of the starting models describing the same underlying 
physics of depletion interactions. 
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Figure 4: Reduced equilibrium density of chain segments X as a function of the 
cohoid volume fraction y. N/n = 1.6 10^ a/D = 15 , Ip/D = 25 , t = 0.39. 



4 Discussion 

The model we propose in this paper allows us to describe in a simple way 
both the depletion interactions induced by condensing agents and the effect of 
correlations between agents on these interactions. Moreover this can be done 
without specifying the kind of depleting agents. Therefore it can be used to 
compare the relative efficiency of spherical colloids and flexible polymers for 
collapsing a single semi-flexible chain. The magnitudes of the induced depletion 
interactions are respectively 



vsa^fj,'{y) 



and 



(18) 



for colloids and polymers. In the case of polymers, the value of the cross virial 
coefficient between a rodlike segment and a polymer coil can be found by requir- 
ing a linear dependence in both the length of the segment (Ip) and the degree 
of polymerization (M).^^ The crossed virial term is then given by 



(19) 



with h the monomer size of flexible polymer. With this expression one can esti- 
mate the strength of the depletion interactions, by using the chemical potential 
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of polymers in a dilute solution: 



^Tfh -yMij,) (20) 

where y '-^ c6'^ is the monomer volume fraction. For colloids the equivalent 
expression reads 



•Across _ V a J /qi \ 

vsa^fi'iy) ^ D/a ^ ' 

where y is the colloid volume fraction. Therefore at equivalent volume fraction 
there is a factor M favoring the PEG compared to the colloid, and as a rule of 
thumb flexible polymers should be much more efficient to collapse a semi-flexible 
chain provided that M >> 1. More precisely, the polymers are more efficient if 

M > ^ JJJ^ (22) 

The enhanced efficiency of polymers compared to colloids for chain condensa- 
tion, at comparable volume fraction, is mainly due to: (i) the smaller transla- 
tional entropy compared to the equivalent fluid of monomers and (ii) the greater 
number of interactions between chain segments and monomers. 

As mentioned in the introduction, several recent papers have dealt with 
DNA condensation by depleting agents. In the case of hard sphere depleting 
agents studied by Sear, the analysis performed in the previous section shows 
that a more systematic description of correlation effects between colloids rules 
out simple unconditional DNA collapse above a critical colloid volume frac- 
tion.^ The conditions required for condensation are strongly restrictive, and 
within the range of validity of the present theory it happens only for colloids 
with diameters smaller than that of DNA. Nevertheless, in a real system, this 
statement can be modified by several factors. First, the polyelectrolyte nature 
of DNA has not been taken into account in our model. As shown in detail by 
Ubbink and Odijk,-^ this factor strongly influences the shape and density of 
the condensate. However, since we are mainly interested in the onset of con- 
densation and reswelling, the electrostatic contribution can be included in an 
effective bare second v irial coefficient.!^ Also in the case of prokaryotic cells 
studied by Odijk et aZ. ,123121 the depleting agents are assumed to be negatively 
charged proteins, i.e., similarly charged as DNA, and this should increase the 
value of the cross virial coefficient Vcross between DNA segments and proteins, 
and therefore increase the strength of the depletion interactions according to 
Eall7l The range of colloid sizes required to condense DNA should be corre- 
spondingly larger. The existence of a limited range of volume fractions that can 
condense DNA should hold a priori; too high a concentration of proteins should 
not be able to condense DNA, a fact that has not been realized so far. Apart 
from the depletion effects of those proteins, it is also no t clea r whether DNA can 
be considered to be in a good solvent in the cytoplasmjEDSI this might enhance 
the depletion interactions between DNA segments induced by the proteins. 
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The ^/'-condensation of DNA has also been described recently by comparing 
the free energy of the system in a coil state and in a condensed state.^ This 
approach neglects completely the nature of the collapse transition and therefore 
depletion interactions, but still describes qualitatively both the collapse and the 
reswelling of DNA in PEG solutions. The correlations between PEGs, incorpo- 
rated through polymer scaling laws in both states of DNA (coil and globule), 
are responsible for the reswelling, as emphasized in our work; the underlying 
physics is the same, i.e., a balance between the loss in translational entropy of 
PEG in the collapsed state of DNA and the depletion effects induced between 
DNA segments. 

The main limitation of our work is the choice of the coupling term Eq^ 
which is the simplest way to couple the behavior of the condensing agents and 
the semi-flexible chain. In particular we assume a linear dependence on the 
chain segment density; but once the chain is collapsed, one expects this density 
to increase, and therefore higher order terms in X might be relevant. However 
these terms are expected to favor reswelling, so that there is no need to include 
them if one is only interested in understanding the qualitative behavior of the 
system and locating the onset of condensation and reswelling. Some of these 
terms -e.g. X'^y - can indeed be thought of as renormalizing the bare second 
virial of chain segments. On the other hand, if higher orders in depleting agent 
concentration are included in the coupling term, the depletion effect is enhanced. 
It is hard to evaluate the precise prefactors of such terms, which influence the 
degree of depletion enhancement. The effect of the cross virial term Xy"^, for 
example, can be qualitatively understood by an increase of the effective size 
of the spheres; a third-order cross virial implying one chain segment and two 
spheres has a very similar effect on the free energy as a second cross virial 
between a chain segment and a bigger sphere. More generally we believe that 
the effect of all higher order coupling terms X"y^ will not change the qualitative 
picture drawn from our model. 

This conclusion is further justified by a recent study on mixtures of hard 
rods and hard spheres based on an integral-equation formalism, which shows 
very similar results for the effective second virial coefficients between rodsj^ 
in the case of rods dissolved in a hard sphere fluid, the effective second virial 
between rods is a non-monotonic function of the hard sphere volume fraction. 
Moreover, the volume fraction ym where the effective second virial coefficients 
between chain segments is minimum, does not depend on the the ratio a/ D, a 
feature that is provided also by our model even though the two models differ 
in the precise estimate of the minimum. Within the framework of integral- 
equation theories, the non-monotonic behavior is explained by the potential of 
mean force, i.e. the depiction potential between two rods that shows a growing 
repulsive barrier to the short-range attractive well associated with the depletion 
interactions, as the sphere volume fraction is increased, similar to the case of 
hard sphere mixtures already mentioned. Moreover, these authors observe, as 
we do, changes of sign of higher order effective virial coefficients that limit the 
practical use of such a virial expansion. Since integral-equation theories take 
into account higher order crossed correlations which we do not, it seems that our 
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choice of the lowest-order coupling term suffices to treat the relevant correlations 
in the system, namely the correlations between condensing agents. Similarly, 
earlier integral-equation studies, also by Schweitzer and coworkers, on the de- 
pletion interactions between proteins immersed in a polymer solution, exhibit 
a non-monotonic dependence of the effect ive pr otein second virial coefficient as 
the concentration of polymers is increased!^^'^ Our present formalism could be 
generalized to compute analytically this effective depletion interaction between 
spherical colloids immersed in a polymer solution. The non-monotonic depen- 
dence of the second virial coefficient can then be explained by the correlations 
in the polymer solution. 

In the numerical analysis presented in section |21 we chose also the simplest 
form for the free energy of the pure semi- flexible chain. This is not really a 
restrictive assumption, because taking into account more correlations between 
chain segments should contribute also to an early reswelling. 

5 Conclusion 

We propose in this paper a simple model for semi-flexible chain condensation by 
neutral depleting agents. It is shown that the kind of agents need not be specified 
to obtain the main features of the induced depletion interactions between chain 
segments. In particular, this allows us to address the role of condensing agent 
correlations on the effective interactions. In the case of polymeric depleting 
agents like PEG, this effect has already been taken into account by Grosberg 
et aZ.p^ the effective interaction between DNA segments is a non-monotonic 
function of the monomer volume fraction. Accordingly, we expect a reswelling 
of the DNA at high enough monomer concentration, and a lower bound for the 
PEG length required to condense DNA. In the case of spherical colloids, our 
model predicts the same feature, unlike recent work where correlations have 
been neglected. The effective interactions are non-monotonic as a function of 
colloid volume fraction, and this again imposes strongly restrictive conditions 
on the size of colloids required to condense a semi-flexible chain. For colloids 
interacting only through their hard cores, no realistic range of colloid size is 
found. Nevertheless, if the quality of the solvent (aqueous solution) for the 
chain is decreased, the range of colloid size should be experimentally accessible. 
We gave in section |21 an example of such a situation by solving the equilibrium 
equations within an SPT description of the colloids. Our model can also describe 
in principle the example of DNA condensation by negatively charged proteins 
considered by Odijk; the only input required is the chemical potential of the pure 
protein solution. In this case, one expects the cross virial term between DNA 
segments and proteins to be larger than the pure hard core value considered 
here, due to electrostatic interactions; the strength of the depletion interaction 
and therefore the condensing range of colloid size are increased. We hope that 
this work will motivate further experimental studies to determine whether DNA 
collapses in the presence of small colloids (very likely proteins). 

The qualitative results obtained in this article are confirmed by recent inte- 
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gral equation treatments of similar systems P^SEZI The analytical model pro- 
posed in this paper provides a simple alternative description of the effect of 
correlations between depleting agents on the interactions mediated by them. 
This approach can be seen as complementary to the integral equations theories, 
since it allows us, in a simple and direct way, to investigate the influence of 
solvent quality and colloid charge, etc., on the phase behave of this system. 
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